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ABSTRACT 

We isolated the anomalous part of the cosmic electron-positron flux within a 
Bayesian likelihood analysis. Using 219 recent cosmic ray spectral data points, 
we inferred the values of selected cosmic ray propagation parameters. In the con- 
text of the propagation model coded in GalProp, we found a significant tension 
between the electron positron related and the rest of the fluxes. Interpreting 
this tension as the presence of an anomalous component in the electron-positron 
related data, we calculated background predictions for PAMELA and Fermi-LAT 
based on the non-electron-positron related fluxes. We found a deviation between 
the data and the predicted background even when uncertainties, including sys- 
tematics, were taken into account. We identified this deviation with the anoma- 
lous electron-positron contribution. We briefly compared this model independent 
signal to some theoretical results predicting such an anomaly. 



Subject headings: astroparticle physics — cosmic rays — Galaxy: general 
diffusion — ISM: general — methods: statistical 



1. Introduction 

Over the last few decades observations of cosmic rays established an increasingly sig- 
nific ant and puzzling d eviation from theoretical pr edictions. Se yeral experiments, such as 
TS dGolden et al.lll994h . AMS (Ulcaraz et al.lboooh. CAPRICE faoezio et al.l lioOlh. MASS 
rtGrimani et alJbooj ). and HEAT faarwick et al.lll997l : feeatty et al.ll2004f ) provided a hint 
of an excess of high energy positrons in our locality. Recent measurements of the PAMELA 
satellite confirmed these suspicions by establishing an excess in the posit ron fraction over 
the theoretical predictions for energies above 10 GeV ( lAdriani et al.l 120091 ). The PAMELA 
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data appear to significantly deviate from the background predic tions even when sizeable ex- 



perimental and theoretical uncertain ties are taken into account (IDelahaye et al.ll2009bl 12010 
MertsdJlioia belahave et aJboilh . 



A possible exces s in the electron-p ositron sum was also indicated by AMS (|Aguilar et al. 



2002|), PPB-BETS flTorii et al.l l2008h . and HESS (lAharonian et all 120081 120091 ) . The ex- 
ceptionally precise measurement of the the local electron+positron flux by the Fermi-LAT 
satelli te, at first glance . seem s to partially confirm the electron+positron excess above 100 
GeV (lAckermann et al.l 120101 ) . The deviation between the Fermi-LAT data (especially the 
2010 release ) and the theoretical backgr ound calculation, produced by the numerical code 
GalProp by iStrong k Moskalenkol (119981 ) appears to be significant. These results were re- 
cently confirmed by the PAMELA collaboration which measured the cosmic ray electron flux 
in a similar energy range and found it to be consistent with the Fermi-LAT data. 

The deviation between the measurements and the predicted backgrounds prompted nu- 
merous attempts to exp lain it by invoking; new physics ranging from modification of the 
cosm i c ray propagation dStawarz et al. 



2009; Hu et al 
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summarizes the present situation of these speculations. 

The existence and statistical severity of the electron-positron anomaly depends on the 
theoretical prediction of the cosmic ray background. While the origin of the cosmic rays is 
not fully understood their local observation, coupled with other astrophysical measurements, 
enables us to build and constrain a model of particle production and propagation in our 
Galaxy Such a model is based on the relatively well understood features of particle diffusion 
within the Milky Way. The diffusion is described by the transport equation, subject to an 
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initial source distribution and boundary conditions. The local electron and positron fluxes 
are calculated by solving this transport equation. Besides the lack of precise knowledge of the 
cosmic ray sources, the background prediction is challenging because the propagation model 
has numerous free parameters, such as the convection velocities, spatial diffusion coefficients, 
or momentum loss rates. 

Motivated by possible new physics buried in the Fermi-LAT data, in this work we 
attempt to determine the size of the anomalous contribution in the cosmic electron-positron 
flux. Our strategy involves the following main steps. 



• Finding the cosmic ray propagation parameters that influence the electron-positron 
flux measured by Fermi-LAT and PAMELA the most. 

• Subjecting cosmic ray data, other than the Fermi-LAT and PAMELA electron-positron 
measurements to a Bayesian likelihood analysis, to determine the 68 % (1-0") credibility 
regions of the relevant propagation parameters. 

• Calculating a background prediction, with uncertainties, for Fermi-LAT and PAMELA, 
based on the determined 1-a credibility regions of the propagation parameters. 

• Subtracting the background prediction from the Fermi-LAT and PAMELA measure- 
ment to isolate the anomalous part of the spectrum. 



Since in the process of the likelihood analysis we determine the uncertainty of the electron- 
positron background, we can also quantify the statistical significance of the deviation between 
the cosmic ray data and the theoretical background calculation. 

When contrasted with the earlier literature, our work contains two main novel results. 

1. Demonstration of a significant tension between the electron-positron related and the 
rest of the cosmic ray data in the context of the propagation model coded in GalProp. 

2. The extraction of the anomalous part of the electron-positron flux . We were able to 
obtain these results because we use more data than other similar studies Maurin et al. 2001 



2002: Maurin et al. 2010: Putze et al. 2010: Lin et al. 2010: Trotta et al. 



20101). We used the 



numer ical code GalProp in the Bayesian framework, extending the analysis of iTrotta et al. 
(120 lOl ) to quanti f y the uncertainty in the background contribution of cosmic ray spectra. 
Unlike iLin et al.l ( 120101 ) we do not use gamma ray data because some components of the 
gamma ray flux are thought to be affected by the same (or similar) anomalous contributions 
as the electron-positron flux. Leaving out the calculation of gamma ray propagation also 
speeds up our numerical calculations. We decided to include gamma ray data in our analysis 
at a later stage. 
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While the numerical analysis by lTrotta et al.l ( 120101 ) is very similar to ours the choice of 
the free diffusion and nuisance parameters are different. More importantly, th e use of sub- 



stanti ally more experimental data (219 data points altogether compared to 76 in lTrotta et al. 



( 120 lOh ) . enables us to constrain the background prediction well enough to isolate the anoma- 
lous part of the e + e~ flux. The experimental data we use, come from multiple instruments, 
over a wide energy range, as discussed in section [331 



2. Cosmic ray propagation 



Cosmic rays are highly energetic particles which have their origins locally and remotely in 



the visible universe ( Ginzb urg fc Syrovatskiil Il964j ; iBlandford fc Eichlerl 119871 : iStawarz et al. 



20101 ; lAharonian et al.ll201ll ). They are divided into two main categories: primary and sec- 
ondary. Primary cosmic rays are particles that are accelerated by astrophysical objects, such 
as supernova remn ants. These cosmic rays interact with interstellar matter to create sec- 
ondary cosmic rays (IBlandford fc Eichlerlll987l ; iDelahaye et al.ll2009bl ; iNakamura et al.ll2010l ; 
Aharonian et al.ll201ll ). The majority of cosmic ray electrons, for example, are likely to origi- 
nate from supernova remnants, while cosmic ray positrons are believed to be mainly produced 



via secondary production processes such as spallation and nucleosynthesis (IBlandford fc Eichler 
19871 : lAdriani et aDbood : belahave et~alll2009a[ INakamura et al.lboiol : lAharonian et al.lboilh . 



Co smic ray propagation through the Galaxy is typically quantified using the dif fusion 
model JGinzburg et al Hl9~9~ol : ISchlickeise^booi IPtuskin et al II2OO6I : Istrong et al.ll2007h . Dif- 
fusion of cosmic rays provides a simple explanation for the highly isotropic distribution of 
high energy charged particles and their noticeable retention in the Galaxy. Diffusion results 
from the particle scattering of cosmic rays on random magnetohydrodynamic (MHD) waves 
and inhomogeneities in the Galactic magnetic field. The random nature of the Galactic mag- 
netic field, causes the trajector ies of the cosmic rays to become jumbled, causing them to un- 



dergo a random walk in spac e ( IGinzburg &: Syrovatskiilll964l : IStrong et al.1120071 ; ICotta et al. 



20 111 ; lAharonian et al.ll201lf ). The energy distribution of cosmic rays are modified by energy 
losses experienced by these particles as they propagate through the Galaxy. Energy losses 
arise due to the interaction of the cosmic rays with the interstellar medium and interstellar 
radiation fields. Re-acceler ation due to interstellar shocks and Galactic winds powered b y 
convection also contribute (IStrong fc Moskalenkdll998l ; IStrong et al.l 120071 : Fan et al.ll2010l ). 
while for heavy and unstable nuclei, fragmentation processes also need to be taken into 
account. 



The diffusion model assumes homogeneous propagation of charged particles within the 
Galactic disk (similar to one of the simplest models of propagation called the leaky box 
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model) but it also takes into account cooling effects. The density tp(f,p,t) (per unit particle 
momentum p) of a particular cosmic ray species at a Galactic radius of r can be calculat ed 



solving the cosmic ray transport equation which has the general form (jStrong et al.ll2007l ) 

= ,(,-, P|1 ) + v.(o„w-i ; fl 

Here q(f,p,t) is the source term which depends on the production mechanism of primary 
and secondary cosmic ray contributions. The spatial diffusion coefficient D xx describes the 
scattering of cosmic ray species through turbulent magnetic fields. This propag ation can be 



isotropic or anisotropic, and can be influenced by the cosmic rays themselves (jStrong et al. 



20071 ). Generally, D xx has the form 



D xx = D 0xx p(^^J , (2) 

where /3 = v/c, and R = pc/eZ is the magnetic rigidity of the particles which describes a 
particle's resistance to deflection by a magnetic field. Here Z is the effective nuclear charge 
of the particle, v is its velocity, p is its momentum, e is its charge, and c is the speed of 
light. The energy of high mo mentum cos mic ray electrons and positrons, for example, can 



be approximated by E ~ eR (lHillaslll984j ). The constant exponent 5 indicates the power law 
dependence of the spatial diffusion coefficient D xx . Different regions of the energy spectra can 
have different 5 values, producing a discontinuity in the derivative of D xx . This artificial kink 



is int roduced so that one can fit the B/C ratio data over all energies ( jStrong fc Moskalenko 



1998|). 



In Eq.([T]), V describes the convection velocity which is a function of Galactic radius r 
and depends on the characteristics of the Galactic winds. The convection velocity is assumed 
to have a cylindrical symmetry and increase linearly with height z from the Galactic plane. 
Apart from transporting particles through the Galaxy, convection also causes the adiabatic 
energy losses (or gains) of cosmic rays due to their interaction with the non-uniform flow of 
gas (Galactic winds) with an inhomogeneous magnetic field. This is represented by the term 
V -V. 

Diffusion in momentum space (diffusive re-acceleration) is described by the coefficient 
D pp . This arises from the scattering of cosmic ray particles on randomly moving MHD waves. 
Diffusion in momentum is related to spatial diffusion via 

D D - (3) 
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Here va is the Alfven speed, the parameter w characterises the level of hydromagnetic tur - 
bulence experienced by the cosmic rays in the interstellar medium ( jSeo fc Ptuskinl Il994j ) . 
This parameter is also known as the ratio of MHD energy density to the magnetic field 
energy density. In the last two terms of Eq.(|T]) the parameter ry is the time-scale of the 
fragmentation loss, and r r is the radioactive decay time-scale. 

Observations of galaxies other than ours suggest that cosmic rays are diffus i ng in a 
cylindrical slab, whose height is dependent on the Galaxy itself ( jDelahaye et al.l l2009al ). 
Consequently, the transport equation (generalised and simplified) is solved in a diffusive 
region shaped as a solid flat cylinder. This cylinder, parametrised by the coordinates (r, <ft, z), 
encloses the Galactic plane with height 2L in the z-direction (z G [—L,L]) and a radius of 
R = 20 kpc in the r direction. The solar system is located at (r, (ft, z)= (8.5 kpc, 0, 0), while 
the boundary conditions imposed on this scenario allow the cosmic ray density to vanish at 
the surface of the flat cylinder and particles may propagate freely outside it and escape. The 
rate of energy loss, b(E) is determined by the photon density, strength of the magnetic field 
and the Thomson scattering cross section associated with the cosmic rays. 

To obtain an explicit analytic solution for a particular cosmic ray species, it is possible 
to solve the simplified version of the transport equation, such as 



dif)(E, r) 
dt 



D XX V 2 *P(E, f) + A [b(E)^(E, r)} + q(E, f) 



(4) 



for electrons (jDelahaye et al.ll2008l ). using a Green's function method (IBaltz fc Edsj oil 19981 ). 
However in most cases, that require a realistic description of the astrophysical environment 
which produce the experimentally observed cosmic ray spectra, an analytical solution is not 
possible. Hence a numerical solution is pursued. 

The numerical Galactic cosmic ray propagation package GalProp calculates the prop- 
agation of relativistic charged particles and their diffuse emission produced during their 
propagation through the Galaxy. GalProp solves the propagation equation numerically for 
Z > 1 nuclei, as well as for electron s and positrons on a two dimensional spatial grid with 
cylindrical symmetry in the Galaxy (jStrong et al.ll2007l ). It also has the capability of solving 
the diffusion equation in three dimensions. GalProp starts with the heaviest primary ele- 
ment defined by the user and the propagated solution is used to compute the source term for 
the secondary products of this element. This process is continued until protons, secondary 
electrons, positrons and anti-protons are produced and a steady state solution is obtained. 
The cosmic ray spectrum is used to compute the gamma rays and energy losses such as 
synchrotron radiation experienced by the cosmic rays. These are computed in conjunction 
with realistic maps of the interstellar gas distributions and radiation fields based on cur- 
rent HI and CO surveys and detailed theoretical calculations of the Galactic magnetic field 
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f lStrong et all 120071 ) 



The input parameter file for GalProp has a number of free parameters which are available 
for the author to define. The main free parameters determine the geometry of the model 
(radius, height of cylinder, and grid spacing), the distribution of cosmic ray sources (which is 
usually chosen to represent an even distribution of supernova remna nts), the primary c osmic 
ray spectral shape and the isotropic composition of the sources (IStrong et al.l 120071 ) . the 
spatial and momen tum diffusion coefficients and their dependence on the particle rigidity 
( jGrasso et al.ll2009l ). These can be classified into a number of subsets: the diffusion of cosmic 
ray, the primary cosmic ray sources and radiative energy losses of these primary cosmic rays. 
The diffusion subset is described by the parameters defined above: 



D 0xx ,S,L,v A ,dV/dz. 



(5) 



The most relevant parameters in the primary cosmic ray source subset are: 



T>e e jynuc 
re.fi i i ref 



nucleus _ nucleus 



7 



(6) 



Here 7 e is the primary source electron injection index. This specifies the steepness of the 



electron injection spectrum, dq(p)/dp oc p 1 
a separate injectio n index for nu c lei de fi 
refer the reader to IStrong et al.l ( 120071 ) . 



below a reference rigidity R e re j- There is also 
a separate injection index for nuclei defined by j nudeus above R™f eus - For further details we 



3. Bayesian inference 

Our aim is to isolate the anomalous part of the Fermi-LAT and PAMELA electron- 
positron fluxes. To do this first we need to know the non-anomalous, standard astrophysical 
background contributing to Fermi and PAMELA. To extract this background and to de- 
termine its uncertainty we use the cosmic ray measurements which appear to be consistent 
with the background estimates. First we determine the values of the Galactic propagation 
parameters most favored by this part of the data and extract the uncertainties of these pa- 
rameters. Then we use these parameter values to calculate the e~ and e + background and 
its uncertainty and compare this to the measurements of Fermi and PAMELA. This way we 
are able to isolate the size and uncertainty of the contribution of the (possible) new source (s) 
in the electron-positron related fluxes. 

To extract the values of the propagation parameters P = {pi, ■■■,Pn} favored by the 
experimental data D = {d±, g?m} we utilize Bayesian inference. In the Bayesian framework 
the probability density of a certain theoretical parameter pi acquiring a given value is given 
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by the marginalized posterior probability distribution 

N 

P{ Pl \D) = / V(P\D) J] d Pj , (7) 

where the integral is carried out over the full range of the parameters. According to Bayes' 
theorem the posterior probability density over the full parameter space is calculated as 

V(P\D) = C(D\P)^. (8) 

Here the likelihood function C(D\P) is the conditional probability density of the theoretical 
predictions for the data with given parameter values P. Data independent information on 
the parameter distribution is folded in via the prior distribution V(P) and the Bayesian 
evidence S(D), for our purposes, acts as a normalization factor. 

The likelihood function, in our case, is calculated as 

M 

C ^ D \ P ) = U^^ ex ^-^(D,P)/2), (9) 



where 



\j(D.P)= l^J^l) . (io) 



The log-likelihood xf contrasts the central value of the i-th data point di with the theoret- 
ical prediction t; L for given parameter values P, in terms of the combined theoretical and 
experimental uncertainty <Xj. 

For parameter estimation the Bayesian evidence only plays the role of an irrelevant 
normalization. Nevertheless, it is useful to calculate S(D) when assessing the validity of the 
hypothesis quantified as the theory parametrized by P. The evidence is easily calculated 
using the normalization of the posterior density 



/ 



N 

P(P\D)Y[d Pj = 1. (11) 



This enables us to recast Bayes' theorem in the integral form 



JV 

£{D)= C{D\P)V{P)Y[d Pj . (12) 

J 3=1 
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Once the posterior distribution is known we can determine the credibility intervals for 
each of the parameters. We define a credibility region 1Z X for parameter p^ by the collection 
of minimal sized parameter regions supporting x % of the total probability: 

x= I V( Pl \D)d Pl . (13) 

In plain terms, a 68 % credibility interval is the minimal parameter region that contains 68% 
of the area under the posterior distribution. This region gives the value of the parameter 
at l-a certainty. Combined credibility regions over multi-dimensions of the parameter space 
can be similarly defined as the minimal region satisfying 

x= V(p i ,p j \D)dp i dp j . (14) 

J Tlx 



3.1. Parameter choice 

The calculation of the posterior probability distributions V{pi\D) requires us to nu- 
merically integrate over the parameter region where the (cumulative) likelihood function is 
non-negligible. The CPU demand to reliably sample the posterior density depends on the 
number of free theoretical parameters N and the speed of the numerical implementation. In 
the case of the diffusion model encoded in GalProp the number of free input parameters is 
around a hundred and for a given set of parameters the code runs for several minutes on 
a single CPU. This makes it unfeasible to attempt the calculation of the posterior without 
simplifications. 

Fortunately, both the number of relevant free parameters and the running time can be 
substantially reduced. To reduce the dimension of the parameter space we tested the robust- 
ness of the electron-positron flux against the variation of nearly all individual parameters 
and found that it is mostly sensitive to the following propagation parameters: 

P = {T*- n »"*™,6 1 ,6 2 ,D 0xx }. (15) 

Here 7 e and ^ nucleus are the primary electron and nucleus injection indices parameterizing 
an injection spectrum without a break, Si and 82 are spatial diffusion coefficients below and 
above a reference rigidity p , and D 0xx determines the normalization of the spatial diffusion 
coefficient. 

We found that the electron-positron spectra are fairly insensitive to the rest of the 
parameters. We also found that the electron-positron spectrum is not only sensitive to the 
power law dependence of the spatial diffusion coefficient D xx , but the presence of a kink 
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therein. So following IStrong et al.l (120041 ) we introduce two coefficients 5i and 62 which 
parametrize this power law dependen ce in Eq.([2]) below a nd above a reference rigidity. We 
fix this reference rigidity to 4 Ge V as IStrong et al.l (120041 ) . 



Our calculations also confirmed the findings of a recent study by lCotta et al.l (1201 ll ) that 
the electron-positron flux is sensitive to the change of the Galactic plane height L. Indeed 
Seo fc Ptuskinl (119941 ) have shown that there is a connection between L and D 0xx : 



Oxx 



2c(l - S)L l ~ s 
3nw5(S + 2) 



(16) 



Thu s, varying the cylinder hei ght amounts to the redefinition of D 0xx as also noticed by 
Ref. |Pi Bernardo et al.l (j2011bl ). In the light of this, we fix L to 4 kpc and use Dq xx as free 
parameter. 

We treat the normalizations of the e~, e + , p/p, B/C, (SC+Ti+V)/Fe and Be-10/Be-9 
fluxes as theoretical nuisances parameters. 



nuisance 



{$°-, $g+, $5/(7, ®(SC+Ti+V)/Fei ®Be~10/Be-9} ■ 



(17) 



They are kept free because the electron-positron flux is either directly or indirectly sensitive 
to these parameters. On the other hand, prior information is available for these parameters 
enabling us to reduce them to the nuisance level. Since GalProp calculates normalizations 
based on local cosmic ray measurements, the results of this calculation can be used as a 
guideline to the central values of the nuisance parameters. The uncertainties of the normal- 
izations can be reliably estimated by an initial scan over the full parameter spaced] 



Varying the parameters listed in Eq. ( lT5i) and ( TTTT) , we confirmed the result of lTrotta et al. 



(120 101 ) that the electron+positron flux of Fermi-LAT can be well reproduced by the theo- 
retical calculation. We also found that by ch anging these parame ters the theory can match 
well the latest PAMEL A electron spectrum (jAdriani et al.l boill ) and the latest PAMELA 
positron fraction data (jAdriani et al.ll2010al ). (We defer the discussion of the quantitative 
details to the results section.) This demonstrates that varying the selected parameters gives 
us enough flexibility to fit all the observed features of the electron-positron spectra. 

While the Galactic propagation of GeV or higher energy cosmic rays is relatively well 
understood, the propagation of a few GeV or lower energy electrons and p ositrons in the 



turbulent, magnetized interstellar medium, remains a formidable challenge (IPrantzos et al. 



1 During our analysis of e ± related or other data we found that the posterior for prefers about 10 
% lower normalization than the value GalProp determines. Since these normalizations form part of our 
parameters in our plots we use the posterior normalizations rather than the GalProp ones. 
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2 1 Oh . Local effects, such as so 



ar modulation and the geomagnetic cutoff, significantly affect 



cosmic rays at lower energies (jPesce- Rollins fc collaboration! 120091 ) . Since solar modulation 



effects, based on the force field model, are built into GalProp, we take these effects into ac- 
count by varying the value of the modulation potential in the code. Following iGast fe Schael 
( 120091 ). we assume a charge-sign dependent modulation, t hat is positively and negatively 
charged cosmic rays are modulated differently by the Sun. IGast fe Schaell ( 120091 ) conclude 
that the effect of this charge dependent modulation on (PAMELA) positrons is substantial. 
They also show that the modulation effect on the p/p ratio is comparable to the statistical 
uncertainties. As described in the next section, we absorb this effect in the systematic un- 
certainties of the p/p data. Heavier nuclei (B, C, Sc, Ti, V, Fe and Be) can carry higher 
positive charges than that of the proton, but their charge to mass ratio is still lower. Since 
the modulation potential is proportional to the charge to mass ratio, the modulation effect 
on heavier nuclei is even milder. Considering that we use the ratio of their fluxes, most of 
the modulation effect cancels since they are positively charged. So the modulation effect on 
heavier nuclei can also be safely absorbed in the systematic uncertainties. 

To be able to compare with experimental data, we set the positron (electr on) modulation 
poten tial in GalProp to <j> + = 442 (2 ) MV. These values were determined by (IGast fe Schael 



20091 ) for PAMELA. lUsoskin et al.l (120111 ) showed that the time dependence of the solar 



modulation potential is not substantial over the period of PAMELA'S data taking, and 
about the same average values can be used for Fermi- L AT . We set the rest of the GalProp 
parameters to the values promoted by IStrong et al.l (120041 ) . 



3.2. Statistical and numerical issues 

In order to extract the most favored values of the propagation parameters we have to 
calculate the posterior distribution P(pj|.D) using suitable Bayesian priors V{pj). Assuming 
no prior knowledge justifies the use of uniform priors. Since we have a previous knowledge 
about the order of magnitude of our parameters we use uniform priors for the propagation 
parameters (rather than for some functions, such as log, of them). For the nuisance param- 
eters prior knowledge is available in the form of a scan over GalProp predictions which are 
based on local measurements of cosmic ray fluxes different from those listed in Table [TJ Thus 
for our nuisance parameters we use normally distributed priors. 



When evaluating Oi for the log likelihood in Eq. ffTO]) . following iTrotta et al.l (120101 ). we 
ignore theoretical uncertainties and combine statistical and systematic experimental uncer- 
tainties in quadrature 



2 _ 2 | 2 

^i ^ i, statistical ' ° i, systematic 
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This can be done for Fermi-LAT and the latest PAMELA e~ flux. Unfortunately, systematic 
uncertainties are not available for the rest of the cosmic ray measurements. When this is the 
case, as an estimate of the systematics, we define <7j as the rescaled statistical uncertainty 

a i = a i, statistical I T i- (19) 

For simplicity, in this study, we use the same scale factor Tj for all data points where system- 



atic u ncertainty is not available. To remain mostly consistent with the work of iTrotta et al. 



(120101 ). we set this common scale factor to a conservative value that they use: Tj = 0.2. We 



checked that our conclusions only mildly depend on this choice. 

We note that systematic errors in the data are not necessarily normally distributed 
point-to-point errors. In fact, typically systematic errors are correlated, such as a systematic 
shift in the energy scale, and could be described by various probability distributions other 
than a Gaussian. Unfortunately, these probability distributions are not provided by even 
those experimental collaborations that indicate a confidence interval for their systematic 
errors. In the lack of this information, we use the simplest ansatz which is adopted by most 
authors in the literature. This estimate of the systematic errors is a simplified approximation 
of a more complicated situation. Nevertheless, for astrophysical data it captures the essence 
of systematic uncertainties. After all, the simplest cosmic ray flux is a falling power law 
spectrum. For this case a systematic shift in the energy scale, for example, can be re- 
interpreted as a systematic normalization shift of the spectrum. Part of this shift is absorbed 
by our normalization nuisance parameters and part of it is approximated as Gaussian error@ 



Due to the simplicity of the posterior density (ITrotta et al.l 120101 ) and its relatively 
low dimensionality we sample the parameter space P and P nu isance according to a simple 
algorithm. We select random model points from the parameter space according to a uniform 
distribution for P and normally distributed for Pnuisance- While this sampling technique 
is less efficient than the Monte Carlo based ones it enables us to trivially parallelize the 
numerical calculation. It also allows us to simply check the robustness of our results against 
the change of certain assumptions such as the prior, the scale factor r or the adequateness 
of the sampling. 

The simplicity of the likelihood function and the high number of data points used in this 
analysis also makes convergence testing relatively simple. To test the validity of our results 
we can evaluate an approximate val ue of the posterior means , variances and the evidence 



adopting the procedure described by iTierney fc Kadand (119861 ). To assure the adequacy of 



the sampling we can simply increase the number of samples of the posterior density until the 



2 We thank the referee of our manuscript to point out this issue. 
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numerically calculated evidence is within 5 % of the one obtained by the Laplace method. 
During this procedure we found that to extract the posterior probabilities presented in this 
paper about one million samples of the posterior density were required over the parameter 
space in Eq.( ll5p and (fT7|) . The gathering of this sample consumed about 2 x 10 5 CPU hours. 



3.3. Experimental data 



We included 219 of the most recent experimental data points in our statistical analysis. 
These contained 114 electron-positron related, and 105 Boron/Carbon, anti-proton/proton, 
(Sc+Ti+V)/Fe and Be-10/Be-9 cosmic ray flux measurements. As a number of experiments 
have energy ranges which overlap, we chose the most recent experimental data points in 
those energy ranges. 



For e + + e we used t he mo st recent data from AMS bvlAguilar et al.l (120021) . Fermi 



LAT by lAckermann et all (120101 1 and HESS by lAharonian et all (12008 



20091 ). The energy 



ranges in which we use each experiment is listed in Table [TJ The AMS experiment re- 
ported an excess in high energy positrons for energies greater than 10 GeV. The Fermi-LAT 
collaboration reported a high precision measurement of the e + + e~ spectrum for energies 
from 7 GeV to 1 TeV using its Large Area Telescope (LAT). This spectrum extended their 
previously publish ed electron positron spectrum over an energy range of 20 GeV to 1 TeV 
(lAbdo et al.l 120091 ) and is flatter than results reported by earlier experiments. HESS's at- 
mospheric Cherenkov telescope (ACT) reported a significant steepening of the electron plus 
photon spectrum above one TeV. 

The PAMELA collab oration measured th e flux of the positron fraction e + /(e + + e~), 
between 1.5 and 100 GeV (lAdriani et al.ll2009l ). They observed that this differential positron 
fraction falls slower than expected for energies above 10 GeV. This behaviour is different from 
that of the background of secondary positrons produced during propagation of cosmic rays in 



the G alaxy. Recently PAMELA released the measurement of the e flux alone (lAdriani et al. 
20 111 ) robustly confirming the e + + e~ spectrum by Fermi-LAT. 



Cosmic ray anti-protons can be used to study the production of primary and secondary 
cosmic rays and their transport throughout the Galaxy. Detailed anti-proton spectra requires 
a large number of measurements over a larger e nergy range, with good statisti cs. Previous 



ballo on borne experiments such as CAPRISE98 (IBoezio et al.ll200ll ) and HEAT (IBeach et al. 



200ll ) detected only a small number of an ti-protons with limited statistics. The PAMELA 
satellite experiment (lAdriani et al.l l2010bl ) provided a comprehensive measurement of the 
anti-proton/proton flux ratio for an energy range of 1-100 GeV. PAMELA'S spectrum follows 
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Table 1: Cosmic ray experiments and their energy ranges over which we have chosen the data 
points for our analysis. We split the data into two groups: electron-positron flux related (first 
five lines in the table), and the rest. We perform two independent Bayesian analyses to show 
the significant tension between the two data sets. 



Measured flux 


Experiment 


Energy 


Number of 






(GeV) 


data points 




AMS fAffuilar et al. 2002) 


0.60 - 0.91 


3 


e + + e" 


Fermi-LAT (Ackcrmann et al. 2010) 


7.05 - 886 


47 




HESS (Aharonian et al. 2008. 2009) 


918 - 3480 


9 


e +/(e+ + e ~) 


PAMELA (Adriani et al. 2010a) 


1.65 - 82.40 


16 




PAMELA (Adriani et al. 2011) 


1.11 - 491.4 


39 


anti-proton/proton 


PAMELA fAdriani et al. 2010b) 


0.28 - 129 


23 




IMP8 (Moskalenko et al. 2002) 


0.03 - 0.11 


7 




ISEE3 (Krombel & Wiedenbeck 1988) 


0.12 - 0.18 


6 


Boron/Carbon 


Lezniak & Webber (1978) 


0.30 - 0.50 


2 




HE AO 3 (Eneelmann et al. 1990) 


0.62 - 0.99 


3 




PAMELA (et al. 2008) 


1.24 - 72.36 


8 




CREAM (Ahn et al. 2008) 


91 - 1433 


3 


(Sc+Ti+V)/Fe 


ACE (Davis et al. 2000) 


0.14 - 35 


20 




SANRIKU (Harevama 1999) 


46 - 460 


6 




Wiedenbeck & Greiner (1980) 


0.003 - 0.029 


3 




Garcia-Munoz et al. (1981) 


0.034 - 0.034 


1 




Wiedenbeck & Greiner (1980) 


0.06 - 0.06 


1 


Be-10/Be-9 


ISOMAX98 Hams et al. (2001) 


0.08 - 0.08 


1 




ACE-CRIS (Davis et al. 2000) 


0.11 - 0.11 


1 




ACE (Yanasak et al. 2001) 


0.13 - 0.13 


1 




AMS-02 (Bureer 2004) 


0.15 - 9.03 


15 
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the same trend as other recent anti-proton/proton ratio measurements. The energy range 
over which we use the PAMELA experiment for the anti-proton/proton ratio is listed in 
Table [U 

In comparison to primary /primary or secondary /secondary cosmic ray ratios, stable 
secondary to primary cosmic ray ratios, such as Boron/Carbon and (Sc+Ti+V)/Fe ratio, 
are the most sensitive to variation in the propagation of cosmic rays in the Galaxy. Their 
sensitivity arises from the fact that primary cosmic rays are generated by the original source 
while secondary cosmic rays are created by t he interaction of their primaries with the inter- 
stellar medium (IChilders fc Duvernoisl 120081 ) . Primary/primary and secondary/secondary 
cosmic ray ratios have a low sensitivity to variation in the propagation parameters as 
the denominator and numerator are produced by similar propagation mechanisms. The 
Boron/Carbon and (Sc+Ti+V)/Fe ratio provide an indication (over different energy ranges) 
of the am ount of interstellar material that primary cosmic rays traverse as a function 
of energy ( IChilders fc Duvernoid 120081 ). The experiments used to define the B/C and 
(Sc+Ti+V)/Fe ratio for our analysis are found in Table [TJ 

In conjunction with stable secondary /primary ratios such as the Boron/Carbon ratio, 
unstable isotope ratios such as Beryllium-10/Beryllium-9 can be used to co nstrain the time 
it takes for cosmic rays to propagate through the Galaxy (jMalininl 12004). In this wor k 
we use Be-1 0/Be-9 data from various experiments, such a s ISOMAX98 (IHams et al.l 120011) . 
ACE-CRIS flDavis et alJl200oh . ACE (TYanasak et al.ll200lh and AMS-02 feurgerll2004h . 



Results 



4.1. Is there a cosmic ray anomaly? 



We begin our results section by investigating the question whether the present cosmic 
ray data can be used to justify the existence of an anomaly in the cosmic electron-positron 
spectrum. Both the reality of an anomaly in the PAMEL A e + /(e + + e~) flux and the absence 
of suc h in the anti-proton flux have been qu e stione d by iKatz et al.l (120091 ) and iKane et al. 
( 120091 ). respectively. Recently iTrotta et al.l ( 120101 ) argued that the Fermi-LAT data can 
be well matched by the diffusion model, as encoded in GalProp, simply by adjusting the 
parameters of the propagation model. Their Fig. 8 clearly shows that the Fermi-LAT 
data agree reasonably we ll with the propaga tion model that was their best fit to 76 cosmic 
ray spectral data points. ITrotta et al.l ( 120101 ) also acknowledge that the "positron fraction, 
shown in the right panel of Fig. 8, does not agree with the PAMELA data (Adriani et al. 
2009), but this was expected since secondary positron production in the general ISM is not 
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capable of producing an abundance that rises with energy". In other words, they conclude 
that PAMELA cannot be fitted by simply adjusting the propagation parameters. We take 
this as an important indication that the cosmic ray anomaly is real and requires a detailed 
investigation rather than the adjustment of the propagation model to explain it. 

The hypothesis that the adjustment of the propagation parameters does not solve the 
cosmic ray anomaly is further supported by the fact that not all cosmic ray data can be 
fitted well w i th a s ingle set of these parameters. It is already evident from Fig. 7 and 8 of 



Trotta et al.l (120101 ) that the best-fit of the propagation parameters to the rest of the cosmic 
ray data does not fit well AMS, Fermi and HESS simultaneously. This is exactly what we 
find in our analysis. Our best fit for all cosmic ray data excluding AMS, Fermi, HESS and 
PAMELA data with electron and/or positron fluxes gives a y 2 per degree of freedom of 
0.34. As a consequence the best fit curves all pass through the estimated systematic error 
bands, shown in gray, in Fig. [TJ When this fit is compared to the AMS, Fermi, HESS and 
PAMELA electron and/or positron flux the x 2 P er degree of freedom we obtain is 24, which 
signals considerable tension bordering exclusion. The converse also holds. By changing the 
propagation parameters, we can find an ideal fit for the electron-positron related fluxes with 
X 2 per degree of freedom of 1.0. But for the rest of the cosmic ray data the same fit results 
in a x 2 P er degree of freedom of 3.1, which is a significant pull for 105 degrees of freedom. 
These discrepancies signal a statistically significant tension between the electron-positron 
measurements and the rest of the comic ray data. 

To further investigate the tension, we divide the cosmic ray data into two groups: 
114 measurements containing observations of electron and/or positron fluxes (AMS, Fermi, 
HESS, and PAMELA) and the rest of 105 data points (anti-proton/proton, Boron/Carbon, 
(Sc+Ti+V)/Fe, Be-10/Be-9). We perform a Bayesian analysis independently on these two 
sets of data extracting the preferred values of the propagation parameters. Remarkably, 
we found that we can obtain information about the electron-positron related propagation 
parameters from the rest of the data because the propagation of the cosmic rays is entangled 
for several reasons. First, certain propagation parameters, most importantly for us Dq xx , 
are species independent. Second, the transport equation includes nuclear fragmentation and 
decay, which directly contributes to the secondary electron-positron fluxes. Third, since their 
energy density is comparable to the interstellar radiation and magnetic fields, various species 
of cosmic rays affect each other dynamics. 

Due to the correlations pointed out above certain parameters of the electron-positron 
propagation are constrained even if no electron-positron related data is used in our fit. 
Unfortunately, the injection indices remain virtually unconstrained. In order to fix those 
parameters we resorted to use a minimal amount of information from the electron-positron 
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Fig. 1. — Best fit curves compared to non-electron-positron related data. The curves 
were calculated using the most probable parameter values inferred from the p/p, B/C, 
(Sc+Ti+V)/Fe and Be data. These most probable values correspond to the peak values 
of the posterior probabilities shown in red in Fig. |2j The best fit curves pass through the 
estimated systematic error bands, shown in gray. 
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Fig. 2. — Marginalized posterior probability distributions of propagation parameters listed 
in Eq. (ll5p . The dashed blue curves show results with likelihood functions containing electron 
and/or positron flux data while the likelihood functions for the solid red curves contain only 
the rest of the comic ray data. Shaded areas show the 68 % credibility regions. A statistically 
significant tension between the electron-positron data and the rest is evident in the three 
lower frames. 
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related fluxes. We decided to use data points from the + e + spectrum because (1) they 
span the widest energy range and the end points of the e~ + e + spectrum, partially due 
to their high uncertainty, appeared to agree with the theoretical predictions even before we 
set out to find the most optimal parameter values; (2) in the low energy region they are 
relatively insensitive for solar modulation effects; and (3) because in the mid energy range 
the e~ + e + theoretical prediction develops an insensitivity to the values of the propagation 
parameters (c.f. the distinct bow-tie shape of the theory uncertainty band). 

With this in mind, we included four e ± related data points in the analysis together with 
the p/p, B/C, (Sc+Ti+V)/Fe and Be data. These were the lowest energy point of AMS, the 
highest energy point of HESS, and the 19.40 GeV and 29.20 GeV data points of Fermi-LAT. 
We have checked that our result are robust against this choice and do not bias the final 
conclusion. 

Fig. [2] clearly shows that the two subsets of cosmic ray data are inconsistent with the 
hypothesis that the cosmic ray propagation model and/or sources implemented in GalProp 
provides a good theoretical description. The five frames display the marginalized posterior 
probability densities of our selected propagation parameters. Dashed blue curves show results 
with likelihood functions containing only electron-positron related flux data (AMS, Fermi, 
HESS, and PAMELA) while the likelihood functions for the solid red curves contain only 
the rest of the comic rays (anti-proton/proton, Boron/Carbon, (Sc+Ti+V)/Fe, Be-10/Be- 
9). Shaded areas show the 68 % credibility regions for the parameters. Table [2] shows 
the numerical values of the best fits and the 68 % credibility ranges for the propagation 
parameters. 

Table 2: Best fit values of the propagation parameters and their 68 % credibility ranges. 
Numerical values are shown for both fits: including the electron-positron related cosmic ray 
data only, and including the rest of the data. 



parameter 


Fit for the 
best fit value 


related data 
68% Cr range 


Fit for the rest of the data 
best fit value 68% Cr range 




2.55 


{2.45, 2.60} 


2.71 


{2.54, 2.92} 


nucleus 


1.60 


{1.51, 1.69} 


2.10 


{1.88, 2.92} 


6i 


0.24 


{0.23, 0.26} 


0.06 


{0.04, 0.08} 


6 2 


0.10 


{0.08, 0.12} 


0.35 


{0.32, 0.39} 


D 0xx [xlO 28 ] 


2.17 


{1.85, 2.19} 


11.49 


{8.86,13.48} 



In the first two frames, showing the posterior densities of the electron and nucleus 
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injection indices 7 e and j nucleus ^ there is a mild but tolerable tension between the electron- 
positron related and the rest of the cosmic ray data. The last three frames, on the other 
hand, indicate statistically significant tension between the e + -e~ and the rest of the data. 
The 68 % credibility regions for the spatial diffusion coefficients 8\ and 62 and that of D 0xx 
fall far away from each other when determined using the two different cosmic ray data sets. 
Although it is not shown, it is easily inferred that not even the 99 % credibility regions 
overlap. It appears that by adjusting the cosmic ray parameters we can, indeed, achieve a 
good fit to either the electron-positron related fluxes or to the rest of the data but not to 
both simultaneously. 

Our interpretation of the tension between the electron-positron fluxes and the rest of 
the cosmic ray data is the following. The measurements of PAMELA and Fermi-LAT are 
affected by new physics which is unaccounted for by the propagation model and/or cosmic 
ray sources i ncluded in our calcul ation. We base this hypothesis partly on the earlier quoted 
statement of iTrotta et al.l (120101 ) that "secondary positron production in the general ISM 
is not capable of producing an abundance that rises with energy". The behavior of the 
PAMELA e + /(e + + e~) data is unexpected based on general theoretical principles and when 
it is fit by adjusting the propagation parameters it leads to a bad fit to the rest of the 
data. An anomaly in PAMELA e + / (e + + e~) is also expected to produce an anomaly in the 
Fermi-LAT e + + e~ and the PAMELA e~ spectra. 



We note that the recently released PAMELA e~ flux flAdriani et al.l 1201 if ), which is 
included among our electron-positron related data, considerably increases this tension. We 
checked that without the inclusion of the PAMELA e~ flux the tension is noticeably milder. 
This, and the e ffect of the extra da ta that we use, probably explains why this tension was 
not detected by ITrotta et al. 



4.2. What is the size of the anomaly? 

We attempt to extract the size of the new physics signal, after arriving to the conclusion 
that new physics is buried in the electron-positron fluxes. Based on our findings our work- 
ing hypothesis is that the new physics is affecting the electron-positron fluxes but hardly 
influences the rest of the cosmic rays. Under this hypothesis the cosmic ray propagation 
parameters can be determined from the unbiased data: anti-proton/proton, Boron/Carbon, 
(Sc+Ti+V)/Fe, Be-10/Be-9. This means that we can use the central values and credibility 
regions of the parameters determined using this data to calculate a background prediction for 
all cosmic ray data including the electron-positron fluxes. Once we quantify the background 
we can subtract it from the electron-positron data to see whether there is a statistically 
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significant signal can be extracted. 

In the first step, we use the central values of the propagation parameters determined 
earlier using p/p, B/C, (Sc+Ti+V)/Fe, Be-10/Be-9 to calculate a central value prediction 
for the PAMELA and Fermi-LAT electron-positron fluxes. Then we use all the scanned 
points in the parameter space lying within the 68 % credibility region of all the five scanned 
parameters to establish a 1-a band around this central value. We will refer to this band as 
the 1-a uncertainty of the background. We overlay this uncertainty band on the Fermi-LAT 
electron+positron and the PAMELA electron and positron fraction fluxes. 

Fig. |3] shows the comparison of the measured electron-positron fluxes and their back- 
grounds. Statistical and systematic uncertainties combined in quadrature are shown for 
Fermi-LAT and PAMELA e~, while (r = 0.2) scaled statistical uncertainties are shown for 
PAMELA e + /(e + + e~) as gray bands. Our background prediction is overlaid as magenta 
bands. The central value and the 1-a uncertainty of the calculated anomaly is displayed 
as green dashed lines and bands. As the first frame shows the Fermi-LAT measurements 
deviate from the predicted background both below 10 GeV and above 100 GeV. 

As we shall discuss later, the low energy deviation might be due to the inadequacies 
of the propagation model, so here we concentrate on the deviation between the background 
and the measurements above 100 GeV. In our interpretation this is a weak but statistically 
significant signal of the presence of new physics in the electron+positron flux. Based on the 
difference between the central values of the data and the background a similar conclusion 
can be drawn from PAMELA. Unfortunately, the sizable uncertainties for the PAMELA 
measurements prevent us to claim a statistically significant deviation]^ 

After having determined the background for the electron-positron fluxes, we can sub- 
tract the background from the measured flux to obtain the size of the new physics signal. We 
determine the central value of the signal by subtracting the central value of the background 
from the central value of the data. The 1-a uncertainty of the signal is the quadratically com- 
bined experimental and background uncertainty. Results for the electron-positron anomaly 
are also shown in Fig. El As expected based on the background predictions a non-vanishing 
anomaly can be established for the Fermi-LAT e + + e~ flux, while no anomaly with statistical 
significance can be claimed for PAMELA due to the large uncertainties. 



3 Recently the Fermi collab oration revealed a very preliminary positron fraction measurement nicely con- 
firming the PAMELA results (|Mitthumsirill201 II ) . Even though the Fermi-LAT makes use of only the Earth's 
magnetic field, it appears to have less systematic uncertainties than that of PAMELA. If the officially pub- 
lished Fermi-LAT measurement will indeed reduce the systematic errors to the level of PAMELA'S statistical 
ones, our background will deviate from it revealing a signal also in the positron fraction. 
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Fig. 3. — Electron-positron fluxes measured by Fermi-LAT and PAMELA (gray bands) with 
the extracted size of the electron-positron anomaly (green bands). Combined statistical and 
systematic uncertainties are shown for Fermi-LAT and PAMELA e~, while (r = 0.2) scaled 
statistical uncertainties are shown for PAMELA e + /(e + + e~). Our background predictions 
(magenta bands) are also overlaid. 
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4.3. What is the source of the anomaly? 



Based on the available evidence we can only speculate about the origin of the discrep- 
ancy between the data and predictions of the cosmic electron-positr on spectra. S ince the 



publication of the first PAMELA positron fraction measurement by lAdriani et al.l (120091 ) 
speculation has been abundant. The first obvious assumption is that some aspect of the 
propagation model used in the present calculat ion is insufficient for the proper description of 



the electron-positron fluxes arriving at Earth (IStawarz et al.l |2010| ; iDonato fc Serpicol 12011 



Arakida fc Kuramatall201ll ; iTawfik fc Salehll2010l ). In this case there exists no anomaly in 
the data. One such plausible effect, which is missed by the 2-dimensional calc ulation in Gal- 
Prop , is the spectral hardening of cosmic rays caused by non- steady sources (jCowsik fc Lee 
19791 ; iPohl &: Espositol Il998l ; iPohl et al.l 120031 ). It would be an int eresting exercise to re- 



peat our analysis using a different calculat ion, such as DRAGON by 



(lEyoli et al. 



2008 



Pi Bernardo et al.l 12010), USINE by Maurin et al. 



Maccione et al.l (12010( 1 



(j201lh . PPPC4DMID 



hv hn,IH „. Mmui n r th, ^ nfhwW ,77\Ui to ^j^Tese possibilities! 



Assuming that the propagation model satisfactorily describes physics over the Galaxy 
the nex t reasonable thing is to suspect loc al effects modifying the electron-positron distri- 
bution (jPesce- Rollins &: collaborationll2009l) . Further suspicion falls on the lack of sources in- 

elude d in the calculation ( Di Bernardo et al. 2011b ; de Vega et al.|[2010 ; Blum 2010 ; Frandsen et al 
20101 ). Possible new sources of cosmic rays to account for the anomaly have been proposed 
in two major categorie s. The first ca tegory is known astrophysical objects with unknown or 



uncertain parameters (jLavalld l201lh . These could 



objects in the Galactic centre, etc. dKawanaka et al 



pe supernova remnants, pulsars, various 



2011 



Ka shivama et al.ll201ll ; iPato et al. 



20101 ; I Yuan et al.l 120 lit iGuo et al.ll2011t iDi Bernardo et al.l 1201 lal ). Finally, more exotic ex- 
planations call for new astronomical and /or particle physics phenomena, such as dark matter 



(Pieri et al. 



2011 



Abidin et al. 



2010 



Josan fc Green 



Cholis fe GoodenoughboiO : 



2010 



Cheng et al. 



2011 



Ko Sz Omura 



Palomares-Ruiz fc Siegal-GaskinslboiO ; 



2010 



Cirelli et al. 



2011b 



Anderson 



2010 



Zah ariias et al. 



2010; 



Yang 



2010; 



Borriello et al. 



201oHKvaell2'oio[ 



2010; iKaiivama fc Okadal 



0; 



2011 



Finkbeiner et al. 



2010 



2011 



2010 



Vincent et al. 2010 



201 ll;lBucklev et al. 



Feldman et al.koidlArina et al.ll20ldlcholisll20ld 



Mohantv et al 



Clind I201CH: llshiwata et al.l 120101: iBarger et al 



2010 



Chen et al 



ogan 



2011 



Bell et al.ll2011b 



Cirelli fe Clind bold: iMasina fc Sannino 



2011 



20101: iKang fc Li 



2011 



Hutsi et al. 



Lineros 201 



Kang et aU2011 



2011 



0; 



Porter et al.l l201lt iHutsi et al.l l201ll ; 



201 



Dugger et al. 



Haba et al. 



Carone et al. 



Zhull201lHBell et alJboilaj ) 



Sanchez fc Holdom 



2011 



Zavala et al 



2011 



Ke et al 



2011 



4 We thank the referee of our manuscript for alerting us to some of the possibilities discussed, and for 
suggesting some of the references, in this paragraph. 
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Possible deviations from the predicted background can occur for energies above 100 
GeV, as electron propagation is limited by energy losses via inverse Compton scattering of 
interstellar dust and CMB light, and synchrotron radiation due to Galactic magnetic field. 
This results in a relatively short lifetime and a rapidly decreasing intensity of the cosmic ray, 
as energy increases. Hence a large fraction of the detected electrons and positrons above 100 
GeV are hypothesise d to come from individual nearby source s that are within a few kilo- 
parsecs of the Earth (IDelahaye et al.ll2009at iGrasso et al.ll2009l ). Random fluctuations in the 
injection spectrum and the spatial distribution of those nearby sources produce significant 
differences in the most energetic part of the observed electron and positron spectrum. This 
can be the indication of new physics either from an astrophysical object (s) or dark matter. 

Whatever the source of the anomaly is, if the size of the anomaly can be isolated then 
the source will have to match that size. Fig. H] compares our extracted signal to a few 
attempts to match the anomaly th at we randomly sele cted from the recent literature. The 
first frame shows the prediction of lAhlers et al.l (120091 ) for unaccounted energetic electrons 
and positrons produced by supernova remnants. The top right frame features contributions 
from additional electron-p ositron prima r y sou rces (nearby pulsars or particle dark matter 
an nihilation) calcul a ted bylGrasso et al.l (120091 ) . The bottom left frame contains predictions 
of iBergstrom et al.l ( 120091 ) for anomalous electron-positron sources from dark matter anni- 



hilatio ns. Similar to this, dark matter annihilation contributions suggested by ICholis et al. 
(120091 ) are shown in the last frame. 



The contributions of various new sources typically come with their own (theoretical) 
uncertainties. In some of the cases this uncertainty is unknown thus it is hard to draw 
any conclusion by comparing these speculations to our isolated signal. In the cases where 
the theoretical uncertainty is known, presently it tends to be large enough to prevent us 
from judging the validity of the given explanation. Nevertheless, based on the present 
amount of information, we can already select a few scenarios that are more favored than 
some others. By adding more data to our analysis it is possible to shrink the uncertainty 
of the signal. Similarly, in most cases, the theoretical model of a given new source can be 
constrained further producing a narrower prediction. With time, more data and more precise 
calculations the various suggestions of the cosmic electron-positron anomaly can be ruled 
out or confirmed. 



5. Conclusions 



We subjected a wide range of cosmic ray observations to a Bayesian likelihood analysis, 
motivated by the possibility of new physics contributing to the measurements of PAMELA 
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Fig. 4. — Comparison of the signal extracted in this work to potential explanations of 
the electron-positron cosmic ray a noma ly . The various theoretic al pr edictions come from 



Ahlers et all (120091 ). iGrasso et all (120091 ). iBergstrom et all (120091 ) and ICholis et all (120091 ). 



Presently the comparison is fairly inconclusive but with more data it will be possible to 
shrink the uncertainty in the determination of the signal. Then various suggestions can be 
confirmed or ruled out. 
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and Fermi-LAT. In the context of the propagation model coded in GalProp, we found a 
significant tension between the electron-positron related data and the rest of the cosmic ray 
fluxes. This tension can be interpreted as the failure of the model to describe all the data 
simultaneously or as the effect of a missing source component. 

Since the PAMELA and Fermi-LAT data are suspected to contain a component unac- 
counted for in GalProp, we extracted the preferred values of the cosmic ray propagation 
parameters from the non-electron-positron related measurements. Based on these parameter 
values we calculated background predictions, with uncertainties, for PAMELA and Fermi- 
LAT. We found a deviation between the PAMELA and Fermi-LAT data and the predicted 
background even when uncertainties, including systematics, were taken into account. Inter- 
preting this as an indication of new physics we subtracted the background from the data 
isolating the size of the anomalous component. 

The signal of new physics in the electron+positron spectrum was found to be non- 
vanishing within the calculated uncertainties. Thus the use of 219 cosmic ray spectral data 
points within the Bayesian framework allowed us to confirm the existence of new physics 
effects in the electron+positron flux in a model independent fashion. Using the statistical 
techniques we were able to extract the size, shape and uncertainty of the anomalous con- 
tribution in the electron+positron cosmic ray spectrum. We briefly compared the extracted 
signal to some theoretical results predicting such an anomaly. 

The authors are indebted to R. Cotta, A. Donea, J. Lazendic-Galloway, Y. Levin, A. 
Mazumdar, N. Sahu and T. Porter for invaluable discussions on various aspects of cosmic ray 
physics and likelihood analysis. KA is thankful to P. Chan for help with issues of parallel 
computing, to M. Jasperse for assistance in Mathematica programming, to J. Lazendic- 
Galloway, T. Porter and to A. Vladimirov for help with GalProp. This research was funded 
in part by the Australian Research Council under Project ID DP0877916 and in part by 
the Project of Knowledge Innovation Program (PKIP) of Chinese Academy of Sciences, 
Grant No. KJCX2.YW.W10. The use of Monash University Sun Grid, a high-performance 
computing facility, is also gratefully acknowledged. 
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